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Abstract 

The present talk summarizes the recently progressed state of a systematic re- 
evaluation of cosmological models that respect the presence of inhomogeneities. Em- 
phasis is given to identifying the basic steps towards an effective (i.e. spatially av- 
eraged) description of structural evolution, also unfolding the various facets of a 
"smoothed-out" cosmology. We shall highlight some results obtained within Newto- 
nian cosmology, discuss expansion laws in general relativity within a covariant fluid 
approach, and put forward some promising directions of future research. 



1 Motivation and Results in Newtonian Cosmology 
1.1 Bridging the Gap 

Does an inhomogeneous model of the Universe evolve on average like a homogeneous solution of Einstein's 
or Newton's laws of gravity ? This question is not new, at least among relativists who think that the 
answer is certainly no, not only in view of the nonlinearity of the theories mentioned |2r| . The problem 
was and still is the notion of averaging whose specification and unambiguous definition turned out to be 
an endeavor of high magnitude, mainly because it is not straightforward to give a unique meaning to 
the averaging of tensors, e.g., a given metric of spacetime. This problem seems to lie in the backyard 
of relativists who, from time to time, add another effort towards a solution of this problem. On the 
other hand, the community of cosmologists "should" locate exactly this problem at the basis of their 
evolutionary models of the Universe. Although there are many exceptions (e.g. ||], 0, Q, 

p7| , p7) |; [ p3t , pi} ], fl44| , a certainly incomplete list), most researchers in this field are drawn back to the 
historical development of cosmologies starting with Friedmann, Einstein and de-Sitter at the beginning 
of the last century. Despite the drastic changes of our picture of structures in the Universe on large 
scales, still, the cosmologist's thinking rests on the hegemony of the so-called "standard model" (i.e. the 
family of FLRW models for homogeneous and isotropic matter distributions). This standard model, up 
to the present state of knowledge, explains (or better is employed to explain) a wide variety of orthogonal 
observations, and it is therefore hard to beat due to its (suggestively) established status of resistance 
against observational tests. Therefore, most discussions in this field are based on the vocabulary of the 
standard model, aiming to constrain its "cosmological parameters" , often on the basis of observations of 
structure in the regional Universe that is very different from homogeneous and isotropic. 

Bridging the gap between an involved mathematical problem of general relativity and the practical 
modeling of cosmological dynamics is possible, if ambiguities of averagers could be removed, and results 
related to contemporary discussions in observational cosmology. In the sequel we shall follow a line of 
thought that will match this need and could intensify work directed to mastering an inhomogeneous 
spacetime. 



1.2 Setting the Pace 

Averaging procedures can be defined in a vast variety of ways (see, e.g., a recent summary by Stoeger et 
al. [^5)). A smoothing-out operator can live on different foliations of spacetime; it can leave an averaged 
field space-dependent, thus operating on a given spatial scale; it can also smooth out all inhomogeneities. 
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It could involve statistical ensemble averaging. It can act on scalars, vectors or tensors; it can smooth 
matter variables or, most elegantly, the geometry of spacetime itself [ plj , p2| . It can average spacetime 
variables and not merely spatial variables, etc. To remove ambiguity here, it is best to work within the 
framework of Newtonian cosmology as a first step of understanding, since there the choice of foliation is 
not a problem and spatial averaging can be simply put into practice by Euclidean volume integration. As 
a second step, we confine ourselves to scalar variables. This allows us to get at least some corresponding 
answer in general relativity, since the averaging of scalars is, for a given foliation, a covariant operation. 
We therefore propose to look at the simplest (mass-conserving) averager as follows. 



1.3 A Newtonian Averager 

Consider any simply-connected spatial domain in Newtonian spacetime. With {-) v we denote spatial 
averaging in Eulerian space, e.g., for a spatial tensor field ^4(x, t) — {^(x, t)} we simply have the 
Euclidean volume integral normalized by the volume of the domain: 

(A) v (t) = ^ f) J^d 3 xA( X ,t) . (1) 

(Here, x are non-rotating Eulerian coordinates.) Since, with this averager, the averaged field is only 
time-dependent, the space-dependence is only implicit by a functional dependence on the domain's 
morphology and position. We shall evolve the domain in time by preserving its mass content. This is a 
natural assumption, if we also want to extend the domain to the whole Universe. 

Now, consider the volume of an Eulerian spatial domain T> at a given time, V — d 3 x, and follow 
the position vectors x = f(X,i) of all fluid elements (indexed by the Lagrangian coordinates X) within 
the domain. Then, the volume elements are deformed according to d 3 a; = Jd 3 X, where J is the Jacobian 
determinant of the transformation from Eulerian to Lagrangian coordinates. The total rate of change of 
the volume of the same collection of fluid elements may then be calculated as follows: 

= Id t / d 3 X J = i / d 3 X d t J=±- [ d 3 X 9 J = (9) v =: 3H V , (2) 

where dt is the total (Lagrangian) time-derivative, 9 — V • v = j the local expansion rate of a given 
velocity model v(x, t); T>i denotes the initial domain found by mapping back V with the help of f _1 
(provided f _1 exists), and Hx> = d-p/a-p naturally defines an effective Hubble-parameter on the domain. 



1.4 Commuting Averaging and Evolution 

From the point of view of the standard model of cosmology it makes no difference, if we smooth the initial 
inhomogeneities (say, at a time in the matter dominated epoch) and then evolve the smoothed data with 
a FLRW solution, or if we evolve these inhomogeneities until present and then smooth the distribution. 
The conjecture is held that both ways produce the same values for the characteristic parameters of a 
FLRW cosmology. 

However, averaging, as defined above, and evolving inhomogeneities are non-commuting operations. 
To see this we have to notice that the total (Lagrangian) time-derivative does not commute with spatial 
averaging in Eulerian space. For an arbitrary tensor field A we can readily derive, with the help of the 
above definitions, the following Commutation Rule ]l8| ]: 

d t (A) v - (d t A) v = (A9) v - (9) v (A) v . (3) 

This tells us that exchanging the operators for averaging and time-evolution produces a source due 
to the presence of inhomogeneities (that we will discuss to consist of positive-definite fluctuations). As 
an example consider the expansion rate itself. Setting A — 9, we get 

d t (9) v - {d t 8) v = (9% - {9% = ((9 - {9) v f) v > , (4) 

i.e., the source is the averaged mean square fluctuation of the expansion rate. It vanishes for the case 
where the local rate equals the global one, which is true for homogeneous-isotropic matter distributions. 
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Figure 1: The evolution of an initial domain T>\ within the Hubble flow of a FLRW cosmology with 
global scale factor a(t) (left) and for a general inhomogeneous setting (right). 



1.5 Formulating the General Expansion Law of Newtonian Cosmology 

As suggested by Eq. (||) we may measure the effective expansion of a portion of the Universe with the 
help of the rate of volume change. Let us introduce an effective dimensionless scale-factor av via the 
domain's volume V(t) = \V\ and the initial volume V\ — V{t\) — \T>\\ (compare Fig. |l|): 

av{t)=( V -P-Y , Le. , H V = ^ . (5) 

For domains T> with constant mass Mp, as for Lagrangian defined domains, the average mass density 
evolves as: 



, (g(*0) Pl M v 

{q) V = — 73 = -sTT- (6) 



Looking at the Commutation Rule, Eq. (Q) , we appreciate that this rule furnishes a purely kinematical 
relation provided we give a velocity model. Dynamics enters by saying how this velocity model is generated 
by gravity. In Eq. (Q) we can express most of the terms already through the scale-factor and its time- 
derivatives. The only unknown is the evolution equation for the local expansion rate. This is furnished 
by Raychaudhuri's equation (which employs the Newtonian field equation V • g = A — 4irGg for the 
gravitational field strength of dust matter g = v). 

Inserting Raychaudhuri's equation, 



9 = A- AirGp - ^9 2 + 2(lu 2 - a 2 ) ; a := y ^o- lj a lj , lu := y ^UijWiij , (7) 

into Eq. (0) , we find that the scale-factor a-p in general obeys the dynamical expansion law : 

3— + 4vrG (g) v - A = Q v , (8) 
av 

with Newton's gravitational constant G, the cosmological constant A, and the source term Qv, which 
we may call "backreaction term" , since it measures the departure from the standard model due to the 
influence of inhomogeneities. 

Qv depends on the kinematical scalars, the expansion rate 9, the rate of shear er, and the rate of 
vorticity oj featuring three positive-definite fluctuation terms: 

Qv = l((° 2 )v-(e)v) + H" 2 )v-H- 2 )v ■ 0) 
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1.6 The Quintessence of Inhomogeneities and Anisotropics 

It is easily verified that Qx> = is a necessary and sufficient condition for having ax> = a(t), with the 
global scale-factor a(t) solving the standard Friedmann equations. Assuming for the moment that the 
averaged Universe would be described by the standard model, then the universal expansion rate would 
be determined by four components: the average mass density (commonly conceived to have baryonic 
and non-baryonic "dark" components), the value of the cosmological constant, and the global "curvature 
parameter" k, which arises as an integration constant by integrating the second-order Friedmann equation 
(Eq. (^) for Qt> = 0) in order to get the global expansion rate 3H(t) = 3|. 

Upon integrating the actual expansion law, Eq. (|J), we instead obtain used a different sign 

convention for Qv)- 

3^ + 3% - SttG (g) v - A = 4" / dt' Q v ^-a 2 v (t') , (10) 

a T> a T> a T> Jti 

where kx> enters as a domain-dependent integration constant. The effective Hubble-parameter Hx> = ^ 
is now determined by the previous (now domain-dependent) components, but also by a fifth component 
due to the presence of inhomogeneities. 

Since there are many unresolved conundrums of the standard model (see, e.g., f27fl , EfJ), cosmologists 
are searching for a fifth parameter that may resolve them (see the talk by Prof. Fujii on "Quintessence" 
in the present Proceedings). The use of the letter Q for the "backreaction" is coincidential, but it may 
possibly be a working coincidence. 



1.7 Cosmic Triangle or Cosmic Quartett ? 

The vocabulary of the standard model may be condensed into the picture of a "cosmic triangle" Q, 
consisting of three cosmological parameters, which are constructed by normalizing the mass density, the 
cosmological constant, and the "curvature parameter" by the square of the Hubble-parameter. 

As suggested by Eq. ( [To|) we are led to define an additional dimensionless "kinematical backreaction 
parameter" through 

in addition to the common cosmological parameters: 



8nG(g) v ~ _ _A_ _ 



However, contrary to the standard model, all f^-parameters are now domain-dependent and transformed 
into fluctuating fields on the domain; for Qx> = the "cosmic triangle" is undistorted and the parameters 
acquire their global standard values. 

Comparing these definitions with Eq. ( |i"o| ) we have 

n£ + fi£ + n£ + n% = i , (is) 

i.e., there are four players in the game (or five, respectively, if we split Vl^ into baryonic matter and 
non-baryonic "dark matter" ) . 

In Friedmann-Lemaitre cosmologies there is by definition no backreaction: fig = 0. In this case a 
universe model with fi TO = 1 conserves the matter parameter. Also, the initial value of the "curvature 
parameter" flfc = remains so during the entire evolution. This changes, if "backreaction" is taken into 
account. £1q itself may act as a "kinematical dark matter", or as a "kinematical cosmological term", 
respectively, depending on the relative strength of shear-, epansion- and vorticity-fluctuations. 



As we have learned |19 and as we shall illustrate below, the parameters corresponding to the "cosmic 
triangle" can experience large changes, even in situations when the "backreaction parameter" is seemingly 
negligible, but non-zero. 
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1.8 Construction Principle for Evolution Models 

In order to obtain quantitative results on the value and impact of the "backreaction parameter" , we have to 
employ evolution models for the inhomogeneities. Their implementation seems not to be straightforward, 
since cosmological evolution models (analytical or N-body simulations) are constructed in such a way that 
they evolve an initial power spectrum of density- and velocity-fluctuations on a given global background 
(a solution of the standard Friedmann equations), implying a vanishing "backreaction" by construction. 
However, the framework of Newtonian cosmology offers the possibility of making use of contemporary 
evolution schemes by assigning sense to a global background within the general expansion law of an 
averaged inhomogeneous model. We shall now elaborate on this fact by demonstrating the validity of the 
following Construction Principle: 

The average expansion of a generic inhomogeneous matter distribution on (topologically) closed Newtonian 
space sections is given by the solution of the standard Friedmann equations. 

This principle, although usually not explicitly stated, lies at the basis of any evolution model in 
Newtonian cosmology. It appears to be very restrictive in light of the general relativistic framework (see 
Subsection 2.5). 

In order to proof this principle let us work with the invariants of the gradient of the velocity fieldQ. 
They are expressible in terms of kinematical scalars and can be written as total diver gen ces of vector 
fields, which has been used and discussed in the context of perturbation solutions ( |pT[ , |25|l): 

I( Viii )=V-v = , H(v w ) = iv • (v(V ■ v) - (v • V)v) = u 2 -a 2 + U 2 . (14) 
The "backreaction term" can be entirely expressed in terms of the first and second invariants: 

Qv = 2(II(v id )}v-l(I(vi, j ))l . (15) 

Now, let us formally assume that there exists a global, homogeneous and isotropic reference model 
with scale-factor a(t). We introduce the following variables with respect to this reference background. 
With the global Hubble-parameter H = a/a we define comoving Eulerian coordinates q := x/a and 
peculiar-velocities u := v — Hx as usual. Using the derivative d qj Ui = dui/dqj with respect to comoving 
coordinates we obtain for the first and second invariants: 

I(v itj )=3H+h(d qjUi ) , ll(vij) = 3H 2 + ^-l(d V]Ul ) + ±11 (d qjUi ) . (16) 
The "backreaction term" remains form-invariant (note: ^d qj Ui 



(17) 



Thus, all terms corresponding to the background flow cancel in this expression, and only inhomogeneities 
contribute to "backreaction" . 



Using Eq. (14) we write the "backreaction term" as a volume-average over divergences. Hence, using 



the theorem of Gauss we obtain: 



n 2 



±f dS-(u(V g -u)-(u-V»-f 




(18) 



with the surface dT> q bounding the comoving domain X> 9 , the surface element dS, and the comoving 
differential operator V 9 . 

From Eq. (|l8| ) we directly obtain Qx> = for a domain with empty boundary, e.g., for toroidal space 
sections, or periodic peculiar velocity fields, respectively. In turn, the assumed reference solution is a 
standard Hubble-flow and may serve as a global background. (For more details see fis)] .) 

This establishes the Construction Principle q.e.d. 

2 A comma denotes partial derivative with respect to Eulerian coordinates d/dxi =,i. 
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1.9 Generalization of the Top— Hat Model 

The Construction Principle provides room for employing the language of standard theories of structure 
formation. As an example we apply the general expansion law to domains within a global Hubble-flow. 

Assuming an Einstein-de-Sitter background for this example we subtract a standard Friedmann 
equation, 3^ + AttGqh = 0, from Eq. (|q) , with the background density qh — H^j, and obtain the 
following differential equation for a-p(i) |l9[ : 

a-p a \ 6 a 



\ a-o a J 2 \ a ) 

with (S) v = {{q) v — Qh)/qh specifying the averaged density contrast 5 in V. For Qn = and {8) v = this 
equation simply states that the time evolution of a domain follows the global expansion, ax>{t) = a(t). 
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For Q v = and 1 + (5) v = the evolution of ax> is still of Friedmann type, but with a mass 

different from the background mass. An important subcase with Qx> = and {5) v =^ is the well- 
known spherical top-hat model [|36| . In Eq. (19|) there are two sources determining the deviations from 
the Friedmann acceleration, the over/under-density and the "backreaction term". This shows that, in 
general, the evolution of a Newtonian portion of the Universe is triggered by an over/under-density and 
velocity-fluctuations. 

1.10 An Averaged Lagrangian Perturbation Scheme 

As a second example we apply the Construction Principle to the widely used and well-tested Lagrangian 
perturbation schemes (see: [^5) for reviews, and Q for a tutorial). 

The trajectory field in the mostly employed "Zel'dovich approximation" j49), |mJ (which can be 
derived as a subcase of the first-order Lagrangian perturbation scheme [ficl ), is given by: 

{ Z (X, t) = a(t) (X + e(i)V lKX)) . (20) 

i/>(X) is the initial displacement field, Vo the gradient with respect to Lagrangian coordinates and £(i) a 
global time-dependent function (given for all background models in ||). 

Given this trajectory field we can already approximate the rate of volume change by the volume 
deformation that is caused by the field. For the effective scale-factor we obtain |p^| : 

(a£ n ) 3 = a 3 (1 + C (Ii)^ + e (IIi) Vi + e . (21) 

A better estimate is to calculate the "backreaction term" from the approximation of the velocity field 
v z = f z and solve the dynamical expansion law, Eq. (||), for av- The "backreaction term" separates 
into its time-evolution given by £(i) and the spatial dependence on the initial displacement field given 
by averages over the invariants of the displacement gradient I,, Hi, and IIIi : 

Qv= 72 x (22) 

(i+^ii^ + e 2 (iii) Ci + e 3 (nii) A ) 
(2 - 2 - + e (e (iiii)^ - 2 - (iii)^) + e (2 (h) Vi (iiii)^ - 1 (iioC 

The numerator of the first term is global and corresponds to the damping factor that also arises in 
the Eulerian linear theory of gravitational instability; in an Einstein-de-Sitter universe £ 2 oc a -1 . The 
denominator of the first term is the volume effect discussed above, whereas the second term in brackets 
features the initial "backreaction" as a leading term and higher-order terms. 

A property of this approximation that renders it very useful as an averaged evolution model should be 
emphasized: it is exact for two orthogonal symmetry assumptions: for the evolution of plane-symmetric 
inhomogeneities, which is a consequence of the known properties of the first-order Lagrangian approxi- 
mation and for the evolution of spherically symmetric inhomogeneities, as will be shown below. 



(i 
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Figure 2: The evolution of the "cosmological parameters" fi^ (solid line), f2f (short dashed), and 
ilq (long dashed) in an expanding (under-dense) domain with initial radius 0.5Mpc (a scaled radius of 
lOOMpc today). The dotted line is marking fi£ + Of + fl^. The upper left plot is for one-c, the upper 
right plot is for three-cr fluctuations (assigned to the expected amplitudes of the initial invariants in a 
Gaussian random field of standard Cold-Dark-Matter fluctuations). The lower plots show the result for 
a collapsing domain. We see here that the "backreaction parameter" is quantitatively more important 
than in an expanding domain - an interpretation of this result will be offered below: a collapsing domain 
experiences more drastic changes in shape. 
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Employing the averaged Lagrangian scheme we can quantify the impact of backreaction on the 
domain-dependent "cosmological parameters" . For details the reader may look at a recent work [ ]19| . 
Here, I would like to mention the remarkable result that for an accelerating, i.e., under-dense region, Uq 
may be numerically negligible as seen in (Fig. |^, upper plots) , but dramatic changes in the other param- 
eters are observed. For a collapsing domain with a present-day radius of lOOMpc the mass parameter of 
the domain may even differ by more than 100% from the global mass parameter (Fig. ^ lower plots). 

1.11 Newton's Iron Spheres 

Looking at the general expansion law, Eq. (||), the careful reader may object that Friedmann's equations 
also hold, if we study the motion of a spherically symmetric domain, although the matter distribution 
inside the sphere may be inhomogeneous. This fact is known as Newton's "Iron Sphere Theorem". Let 
us show that the general expansion law respects this theorem. 

We note that, for a spherically symmetric distribution of matter within the domain, the invariants of 
the velocity gradient, averaged over a ball Br, obey relations resulting in p9| : 

<II(«ij)> Bj| = I (I(«ij)>L , and (III(v i , j )) BR = ± (L(v id ))% B . (23) 

Inserting this into the "backreaction term" , Eq. (|l^) , shows that Q g pherlcal = in accordance with 
Newton's theorem. The inhomogeneous model discussed in the last subsection also has this property: 
inserting the above expressions into Eq. ( p2| ) (using the proportionality of displacement gradient and 
initial velocity gradient), we also obtain Qg = = Qg' enca . 



1.12 Averaging and the Evolution of Form 

The expansion law, Eq. (||), is built on the rate of change of a simple morphological quantity, the 
volume content of a domain. Although functionally it depends on other morphological characteristics of 
a domain, it does not explicitly provide information on their evolution. An evolution equation for the 
"backreaction term" is missing. This fact touches on the problem of closing the hierarchy of dynamical 
evolution equations considered as a set of coupled ordinary differential equations in Lagrangian space. The 
problem of closing such a hierarchy of equations is often considered in the literature and various closure 
conditions are formulated (e.g., jjoj), one of them being the "Silent Universe Model" in general relativity, 
which assumes a vanishing magnetic part of the Weyl tensor 0. Averaging such a hierarchy would 
result in evolution equations for the "backreaction term" and would, with some local closure condition, 
also close the system of averaged equations. Here, we will not pursue this problem further, but instead 
begin to develop the morphological point of view that eventually implies an alternative proposition of a 
morphological closure condition. 

Let us focus our attention on the boundary of the spatial domain T>. A priori, the location of this 
boundary in space enjoys some freedom which we may constrain by saying that the boundary coincides 
with a velocity front of the fluid (hereby restricting attention to irrotational flows). This way we employ 
the Legendrian point of view of velocity fronts that is dual to the Lagrangian one of fluid trajectories. 
Let S(x, y, z, t) = s(t) define a velocity front at Newtonian time t, v = VS. Below, we shall need that the 
three principal scalar invariants of the velocity gradient v^j —: £i 7 - can be transformed into divergences 
of vector fields as written explicitly in Eqs. ( jl4 

Defining the unit normal vector n on the front, n = ±y^§y (the sign depends on the direction of its 
motion), the average expansion rate can be written as a flux integral using Gauss' theorem: 

(e) v =±^d 3 xV-v = ±J^dS-v, (24) 

and, with the surface clement da, dS = nda, we obtain the intuitive result that the average expansion 
rate is related to another morphological quantity of the domain, the total area of the enclosing surface: 

{0)v = Tr\ M^S\ . (25) 



26). 



Inserting the velocity potential also into the other invariants, Eq. (Q) and 

ni(« ili ) = iv-Qv-(v(V-v)-(v.V)v)v-(v(V-v)-(v.V)v)-VvJ , (26) 
and performing the spatial average, we obtain [p"7| : 

(11)^ = 1 f d 3 xll= [ da\VS\ 2 U- (m) D = ~ f d 3 xlll = [ da\VS\ 3 K , (27) 

V Jx> JdV * Jv JdV 

where H is the mean curvature and K the Gaussian curvature of the 2— surface bounding the domain. 
|V5| = 4? equals 1, if the instrinsic arc-length s of the trajectories is used instead of the extrinsic Newto- 
nian time t. The averaged invariants comprise, together with the volume a complete set of morphological 
characteristics related to the Minkowski junctionals of a body: 



W (s) := / d 3 x = V ; Wi(s) :=\ j da ; W 2 (s) :=\ f da H ; W 3 (s) :=\ j da K = % X . 

JV 3 JdV 6 JdV 6 JdV 6 

(28) 

The Euler-characteristic x determines the topology of the domain and is assumed to be an integral of 
motion (x = 1), if the domain should remain simply-connected (a morphological closure condition). 

Thus, we have gained a morphological interpretation of the "backreaction term" : it can be entirely 
expressed through three of the four Minkowski functionals: 

°" (s) ^(t-f ) ■ (29) 

The W Q ; a = 0,1,2, 3 have been introduced as "Minkowski functionals" into cosmology by Mecke 
et al. (p5[) in order to statistically assess morphological properties of cosmic structure. Minkowski 
functionals proved to be useful tools to also incorporate information from higher-order correlations, e.g., 
in the distribution of galaxies, galaxy clusters, density fields or cosmic microwave background temperature 
maps ( j34|, 1 33 , (5J, ^l}; see the review by Kerscher ^32) and ref. therein). Related to the morphology 



of individual domains is the study of building blocks of large-scale cosmic structure |3j| , J42J . 
For a ball with radius R we have for the Minkowski functionals: 

:= f R 3 ; <*(*) := f R 2 ; : = ^-r ; wf( s ) ■■= f • (so) 



Inserting these expressions into the "backreaction term", Eq. (29), shows that (s) = 0, and we have 
confirmed Newton's "Iron Sphere Theorem" once more. Moreover, we can understand now that the 
"backreaction term" encodes the deviations of the domain's morphology from that of a ball, a fact which 



we shall illustrate now with the help of Steiner's formula of integral geometry (see also 35 ). 

Let dao be the surface element on the unit sphere, then (according to the Gaussian map) da = 
RiR 2 daQ is the surface element of a 2— surface with radii of curvature R\ and R 2 . Moving the surface a 
distance e along its normal we get for the surface element of the parallel velocity front: 

da e = (JZi + e)(R 2 + e)da = + ^ ' + R ^ + - da = (1 + e2U + e 2 K)da , (31) 

R\R 2 

where 

11=1(1 + 1-) , , (32) 

2 \R± i?2 / R\R 2 

are the mean curvature and Gaussian curvature of the front as before. 

Integrating Eq. ( |3l| ) over the whole front we arrive at a relation between the total surface area A 

of the front and A s of its parallel front. The gain in volume may then be expressed by an integral of 

the resulting relation with respect to e (which is known as Steiner's formula defining the Minkowski 

functionals of a (convex) body in three spatial dimensions) : 

V e = V + f de'A £ , = V + eA + e 2 f daH + e 3 f da K . (33) 

J0 JdV JdV 
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2 Average Models in General Relativity 



2.1 An Averager for Scalar s 

The key to reducing ambiguity of choosing an averager in general relativity is to confine ourselves to 
scalars, since spatially averaging a scalar field ^ is a covariant operation. However, we have to make 
sure that we understand the word "spatial" physically. Geometrically, we shall introduce a foliation of 
spacetime, but we shall encounter ambiguity in choosing it. Suppose for the moment that we insist on 
spatial averaging and consider for all what follows a foliation of spacetime as given. The simplest averager 
is then immediately written down as follows: 

(^>(X i ,t)) v :=y J Jd 3 X V(X\t) ■ V = J^Jd 3 X , (34) 

with J :— -y/det (<7ij ) , where gij is the metric of the spatial hypersurfaces, and X 1 are coordinates in the 
hypersurfaces, conveniently chosen such that they are constant along flow lines. 



2.2 Foliating Spacetime: Covariant Fluid Approach 

We wish to conserve mass inside a spatial domain of spacetime. Therefore, let us first consider the 
(conserved) restmass flux vector^ 

M M := gu^ ; =0 ; g > , (35) 

where g is the restmass density and the flow lines are integral curves of the 4-velocity u^. We shall 
confine ourselves to irrotational perfect fluids with energy density e and pressure p, which allows us to 
simplify the splitting of spacetime, e.g. the spatial hypersurfaces can be chosen flow-orthogonal in this 
case, the unit normal on the hypersurfaces coincides with the 4-velocity and the shift vector vanishes. 
Irrotationality guarantees the existence of a scalar function S, such that 

■ ~ S ' S (36) 



h 

In general, we identify the magnitude h with the "injection energy per fluid element and unit restmass" , 

h := i±£ , (37) 
Q 

which is related to the relativistic enthalpy r\ := by h = r)/m with m the unit restmass of a fluid 
element, and n the baryon density. Note that de — hdg. For a barotropic fluid we can easily see that e 
is a function of the restmass density only and, hence, ft, is a function of g. h is identical to 1 in the case 
of dust. The magnitude h normalizes the 4-gradient d^S so that u^u^ = —1, 



h = y/-d a Sd a S = u^d^S = S > . (38) 

The overdot stands for the material derivative operator along the flow lines of any tensor field T 
as defined covariantly by T := u^T It reduces in the case of a scalar function in a flow-orthogonal 
slicing (which we want to envisage) to the partial time-derivative multiplied by jt , where N is the lapse 
function (for more details see: |lq| ). 

It can be shown that S is spatially homogeneous. Thus, S(t) and h(X l ,t) play the role of "phase" 
and "amplitude" of the fluid's wave fronts. Foliating the spacetime into hypersurfaces S(t) = const, is a 
gauge condition that is naturally adapted to the fluid itself and, therefore, this foliation can be given a 
covariant meaning (§, § and ref. therein). 

3 Greek indices run through 0...3, while latin indices run through 1...3 as before; summation over repeated indices is 
understood. A semicolon will denote covariant derivative with respect to the 4— metric with signature (— , +, +, +); the units 
are such that c = 1. 
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Figure 3: The evolution of an initial domain T> within a simply-connected tube of spacetime. The 
figure visualizes the restrictions implied by the topological requirement of simple connectivity throughout 
the evolution of the domain. Note that the boundary of the domain may break in the coarse of structure 
formation (occurence of a Legendrian singularity), which may induce topology changes. 



2.3 A Commutation Rule in General Relativity 

The rate of change of the volume V(t) in the hypersurfaces S(t) = const, is evaluated by taking the 
partial time-derivative of the volume and dividing by the volume. Since dt and d 3 X commute (but not 



' and d 3 X !) we obtain: 



-^- = — d 3 Xd t J =Tj d 3 XNJ = — I d 3 XN6J = (N9) v . (39) 



dr ■ N 

~v V J v l " vj v - V JV 

Introducing the scaled (t-)expansion 8 := NO we define an effective (t-)Hubble function in the hypersur- 
faces by 

(9)v = d jY =?> ^ = : 3 H V . (40) 
V av 

With the help of these definitions we readily derive the Commutation Rule: 

8 t - (d t *) v = - (0) v . (41) 

The lapse-weighted quantities have to be introduced only in the case of non-vanishing pressure, since the 
pressure gradient induces deviations from a geodesic flow implying an inhomogeneous lapse function. In 
the much simpler case of a dust matter model, the lapse function is homogeneous and can be chosen equal 
to 1, and the covariant fluid gauge is identical to the comoving and synchronous gauge, which makes the 
correspondence to the Newtonian investigation most transparent. 

In the sequel we shall, for simplicity, discuss the case of a pressure-less fluid only and resume the 
discussion of the more general case thereafter. 



2.4 Expansion Law of General Relativity: Dust Models 

Employing the averaging procedure outlined above and following the line of thought of the Newtonian 
investigation in Section 1, we derive the general expansion law for dust matter in the synchronous and 
comoving gauge. The result will be covariant with respect to this foliation. 

I here give the result derived in jis): the spatially averaged equations for the scale-factor op, respecting 
mass conservation, read: 
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averaged Raychaudhuri equation: 



3^ + 4nG^ - A = Q v ; (42) 
av Via?, 



averaged Hamiltonian constraint: 

'a v \ 2 8nG M v (K) v A _ Q v 



av J 3 ViOjj 6 3 6 



(43) 



where the mass Mt>, the averaged spatial Ricci scalar (JZ) V and the "backreaction term" Qd are domain- 
dependent and, except the mass, time-dependent functions. The backreaction source term is given by 

Q D :=2<II> D -!<I& = !((0- {9) v f) v -2{^) v . (44) 

Here, I and II denote the invariants of the extrinsic curvature tensor that correspond to the kincmatical 
invariants we employed earlier. The same expression (except for the vorticity) as in Eq. (^) follows by 
introducing the split of the extrinsic curvature into the kinematical variables shear and expansion (second 
equality above). 

We appreciate an intimate correspondence of the GR equations with their Newtonian counterparts 
(Eq. (||) and Eq. (|To|)). The first equation is formally identical to the Newtonian one, while the second 
delivers an additional relation between the averaged curvature and the "backreaction term" that has 
no Newtonian analogue. This implies an important difference that becomes manifest by looking at the 
time-derivative of Eq. d43| ) . The integrability condition that this time-derivative agrees with Eq. ([|2]) is 
non-trivial in the GR context and reads: 

d t Q v + 6— Qv + d t (K)v+2—(K)v = ■ (45) 

The correspondence between the Newtonian fc-p-parameter and the averaged spatial Ricci curvature is 
more involved in the presence of a "backreaction term" : 

% - \a 2 v f dt' Q v ^-al(t') = \ ((K) v + Q v ) . (46) 
a x> Ju dt' 6 

The time-derivative of Eq. ( |46| ) is equivalent to the integrability condition Eq. (f45|). Eq. ( fl5|) shows that 
averaged curvature and "backreaction term" are directly coupled unlike in the Newtonian case, where 
the domain-dependent /cp-parameter is fixed by the initial conditions. 

For dust models in general relativity we may therefore introduce dimensionless average characteristics 
that are slightly different from those in Newtonian cosmology, the difference being that the "backreaction 
parameter" is now directly expressible in terms of the "backreaction source term" , and not just as an 
integral expression: 

v _ 8ttGM v v _ A v _ (K) v v _ Q-p 

which also obey Eq. (|l3|). 

The evolution of these parameters is intimately related, unlike the situation in the standard model, 
as may be illustrated by the following equation [15[: 



^+6/f c o2(i-nr-^Q)+^?+2 J ff^(i-^-ng)-3if,,(i-n^-o^-og)(^+n2) = o. (48) 

It is, e.g., a good excercise to show with the help of this equation that an inhomogeneous model (including 
"backreaction"), in which the matter parameter stays 1 like in the standard Einstein-de Sitter cosmos, 
does not exist. 
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2.5 Liberation from Strict Meter 



Let us now hold in for a moment and sort out what the Newtonian and the relativistic expansion laws 
for dust matter distinguish. Their close correspondence bears the temptation of overlooking a crucial 
conceptual challenge for the modeling of inhomogeneities in general relativity. 

Given the Newtonian average model, Eq. (§h, and its quantitative consequences (elaborated in detail 
in fil| ) , we can draw the conclusion that we only have to consider spatial scales that are large enough to 
have a negligible influence from the "backreaction term" . This term, which brought the higher voltage of 
having mastered a generic inhomogeneous Newtonian cosmology, shows no global relevance, and it seems 
that we are drawn back to the previous state of low visibility of the standard cosmological models. 

The key-reason for this outcome is the validity of the Construction Principle. We force a globally 
vanishing "backreaction" by the way we construct inhomogeneous evolution models. (Remember that 
the majority of cosmological N-body simulations and analytical approximations are Newtonian and as- 
sume periodic boundary conditions for inhomogeneities on a FLRW background.) The apparent global 
irrelevance of "backreaction" should not be a surprise, since our evaluations have taken place within a 
periodic cube into which the universe model was forced: the scale of this cube introduces a "strict meter" 
governed by the scale-factor of a standard FLRW universe, the action of fluctuations of inhomogeneities 
is confined to the interior of this cube. 

The relativistic average model not only has a somewhat richer tone by linking the averaged spatial 
Ricci curvature to the "backreaction term" as a result of the Hamiltonian constraint, it also places an 
equal stress on both of them in the following sense. Suppose we let the self-gravitating fluid evolve 
freely and link its structure to the curvature of space sections through Einstein's equations. Let now the 
inhomogeneities evolve out of an almost homogeneous and isotropic state on almost flat space sections. 
Naturally, the inhomogeneities grow in the coarse of structure formation. The average Ricci curvature 
of the space sections, however, is dictated by the fluid structure unlike in the standard model which 
would suggest an evolution inversely proportional to the "strict meter" a(t). In contrast, the value of 
the averaged curvature is an outcome rather than a model assumption. Thus, a genuine property of the 
relativistic model, also globally, is a change of the average curvature that does not follow a predesigned 
global law; a priori it is not constrained on some large scale, because the link to the "backreaction term" 
makes sure that, at any time, the "interior" of a universe model (below this scale) also determines the 
global curvature. In other words: there is no analogue of the Construction Principle. This comment is 
not entirely correct, since further study may reveal that it could be possible to formulate an analoguous 
principle on the basis of topological constraints that may be imposed on the space sections. In any 
case, research should not lack commitment to the challenge of constructing inhomogeneous evolution 
models on curved space sections, since this is evidently the generic case. We note that existing work 
on general relativistic evolution models and their averages relies in most cases on assumptions that are 
closely designed after the Newtonian case in order to rescue standard procedures like periodic boundary 
conditions, decomposition of inhomogeneities into plane waves and set-ups of initial conditions on locally 
flat space sections (compare |37j and comments in |L5| and Jl^|). 

2.6 What Einstein Wanted 

In the previous subsection we advocated a viewpoint that concentrates on the physics of fluids, which 
(actively) determines the geometry of spacetime, in particular its average properties. As a showcase we 
learned in the Newtonian framework that the morphological properties of spatial domains are determined 
by the averages over fluctuations in the kincmatical variables (Subsection 1.12). A similar view may also 
apply to the relativistic context and, most interestingly, also to the global morphology of space sections. 

Turning this around, we may alternatively take the (passive) viewpoint of understanding the global 
world structure as given in terms of geometrical (and topological) conjectures (which may themselves 
be subjected to observational falsification; an example is the possibility of falsifying possible topological 
spaceforms on the basis of microwave background observations (|24|]; see also [Q and ref. therein). 

Einstein's vision of a globally static universe model stands out as a famous example. His "prejudice" 
that the Universe as a whole should be static and all evolution (assigning sense to the notion of time) 
should take place in the "interior" of this world model entails a strong determination of the model's average 
properties. Only his invention of the cosmological term made this vision possible resulting in the beautiful 
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structure of a spherical space in which, apart from topological degeneracies, many characteristics (like 
the radius of the universe model) are practically fixed. Hence, according to what we have learned about 
average models, the fluid's fluctuations are "slaved" to the global predetermined structure. This global 
model even guarantees strong evolution, because fluctuations evolve exponentially in a static cosmos. 

The problem that his model is unstable (within the class of FLRW models) applies also to the other 
expanding candidates, if the class is widened to include inhomogeneities. Einstein's example so serves as a 
guide to think about the average model we have discussed in a geometrical way. The assumption of a static 
model (the general expansion law allows for a static cosmos even without the cosmological term) should 
be taken as an example for the possible slaving of fluctuations to global structural assumptions. Of course, 
other such assumptions are possible. This illustrates how global geometrical and topological constraints 
on spaceforms could provide "boundary conditions" for the evolution of fluctuations. Since, in Newtonian 
cosmology, we are always working on a toroidal space without further questioning this assumption, this 
point of view is not exotic, but plants the seed for possibly fruitful research directions. That these 
directions are mathematically very involved can be made obvious by pointing out the relation between the 
average curvature and its compatibility with topological spaceforms (a spaceform with globally negative 
spatial curvature is not compatible with simply-connected space sections) . Due to the process of structure 
formation the boundary of the domain may break, or it may self-intersect and split into two domains 
(these events are topologically classified in the framework of Legendrian singularities, |l]] and ref. therein). 
Since the spatial metric is linked to the fluid, a singular dynamics could induce topology changes. Contrary 
to the existence of such natural metamorphoses, the averager we have defined relies on the assumption 
that the domain remains simply-connected (Fig. [5]). 



2.7 Expansion Law of General Relativity: Perfect Fluid Models 

As an outlook I briefly give one of the results obtained by averaging a perfect fluid cosmology on hyper- 
surfaces S(t) = const, as explained above. 

The averaged equations can be summarized in the following way |l6f| : 

Let us define effective densities as sources of an expansion law, 



Qv Qv Vt> , 



with the scaled matter sources e :— N 2 e and p := N 2 p. Then, the averaged equations can be cast into a 
form similar to the standard Fricdmann equations: 



(IV 

7-2 



4^G(e cff + 3p off ) = ; (50) 



6H V + {11} V - 167rGe cff = , (51) 



and the integrability condition of Eq. (|5C| ) to yield Eq. ( |5l| ) has the form of a balance equation between 
the effective sources and the averaged spatial (t-)Ricci scalar 1Z := N 2 1Z: 

d t e cS + 3H v {£ctt+Pctt) = j^(d t {K) v + 2H v (n)v) • (52) 

The effective densities obey a conservation law, if the domains' curvature evolves like in a "small" 
FLRW cosmology, {TZ)v = 0, or {71) v K a^ 2 , respectively. 

The expressions Qv and P-p are given in As can be seen above, Qv (the "kinematical backre- 
action" ) acts like a fluid component with "stiff" equation of state, similar to the action of a minimially 
coupled scalar field, while Vv (the "dynamical backreaction" ) is due to the non-vanishing pressure gra- 
dient in the hypersurfaces. 

It is interesting that the "kinematical backreaction" of the inhomogeneities could play the role of a 
free scalar field component that is commonly introduced when early evolutionary stages of the Universe 
are studied. It is also interesting that the expansion law in the case of a minimally coupled scalar field 
source gets particularly simple. This case will be investigated in a forthcoming work p0|. 
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